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Abstract 
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1 Introduction 



It is well known that integrable models have the following factorization property: any 
multi-particle scattering amphtude can be reduced to the product of 2-particle scattering 
amplitudes. The independence of the order in which the successive 2-particle scatterings 
take place leads to the Zamolodchikovs so-called triangle equations for the S-matrix P]. 
In statistical physics, in the context of integrable lattice systems, the same property was 
analyzed earlier by C.N. Yang [2] and R.J. Baxter ^ using the so-called i?-matrices and 
the equations encounted were denoted Yang-Baxter equations (YBE) by L.Faddeev |4j. 

/^From the i?-matrices one can construct the transfer matrices of the integrable lattice 
systems and the Yang-Baxter equations impose sufficient conditions on the i?-matrices to 
ensure that transfer matrices t{u) and t{v) with different so-called spectral parameters u 
and V commute. 

In the analysis of the Zamolodchikovs the factorization property was a consequence of 
an assumed Lorentz invariance in the model. Thus it is natural to ask if it is possible to 
construct an integrable lattice model with a commuting family of transfer matrices t{u), 
where the transfer matrices are not products of only 2-particle /^-matrices, as is the case 
for the presently known integrable models. An obvious, interesting equation is then if 
such lattice models can be associated with a Lorentz invariant continuum theory. 

In this article we will address this problem and will analyze the possibility of con- 
struction of an integrable model with 3-particle i?-matrices which can not be reduced to 
products of three 2-particle i?-matrices. In a previous article [S] we have considered a 
model which appeared naturally as a simplification of the so-called sign-factor representa- 
tion of the three-dimensional Ising model (3DIM) on a dual body centered cubic (DBBC) 
lattice ^. A central point in the construction is the appearance of a two-dimensional ran- 
dom Manhattan lattice (ML). Such a Manhattan lattice appears also in the sign-factor 
model (SFM) of the 3DIM formulated on a cubic lattice [B^, where the hopping of fermions 
from site to site is allowed only along the ML arrows. This directed hopping model can be 
described as a product of 2-particle i?-matrices. However, in the case where the surfaces 
were embedded in a DBCC lattice we discovered the presence of 3-particle i?-matrices 
which can not be reduced to the product of two particle i?-matrices. Graphically they 
could be represented as honeycombs [U^, see Fig lb. 

In this article we go further and construct an integrable model where the transfer ma- 
trix is a product of 2-particle i?2-niatrices and (non-reducible) 3-particle i?3-matrices. The 
model is defined in Section 2. In Section 3 we show that the corresponding Yang-Baxter 
equations ensure that transfer matrices with different spectral parameters commute and 
we present a nontrivial solution. We also obtain an explicit representation of the transfer 
matrix as the normal ordered exponential of a (non-local) Hamiltonian. 

"^This representation of the Ising model is in many respects the natural generalization of the repre- 
sentation of the two-dimensional Ising model as a sum over random walks. Here it becomes a sum over 
certain random surfaces, which can be related to fermionic strings. 
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2 Formulation of the model and definition of the /^-matrices 

In in] the SFM for the 3DIM was defined as a theory of fermions interacting with an 
induced Z{2) gauge field on 2d random surfaces embedded in a three-dimensional regular 
cubic lattice. Fermions were hopping along the directed links of a certain random ML. In 
inj we generalized this construction of a SFM to the case of a 3DIM on a DBCC lattice and 
encountered, as mentioned above, not only the ordinary 2-particle R-matrices but also 
3-particle R-matrices which could not be factorized into the product of three 2-particle 
R-matrices. We initiated the analysis of the SFM by considering the restriction of the 
model to a certain regular ML lattice and investigated it's integrability. 

In this article we study in the same spirit a model on a more complicated regular ML. 
Also this model has 3-particle R^- and ordinary 2-particle i?2-matrices, but due to the 
different structure of the ML their appearance in the transfer matrix is different and we 
will show that the corresponding YBE's, unlike the case considered in jS], have spectral 
parameter dependent solutions. 

The transfer matrix of the model is constructed as the product of two monodromy 
matrices Ti and T2: 

T = trT,T2, (1) 
where the matrices Tj are defined as shown in Fig. 2 in terms of the /^-matrices 

T^ = l[Rsin)R2in)R2{n), (2) 

n 

T2 = l[R2{n)R2{n)R3{n). 

n 

The T-matrices are acting on isomorphic quantum spaces which are product of two- 
dimensional spin-spaces at the horizontal sites l(n), 2(n), 3(n), 4(n),... in Fig. 2, and they 
also act on isomorphic auxiliary spaces constructed from the left and right boundary sites 
in Fig. 2. The trace in is taken with respect to the auxiliary spaces. 

Following the technique of J3| one can express graphically the R2- and i?3-matrices as 
shown in Fig.l and the transfer matrix r as shown in Fig. 2. Alternatively this model can 
be formulated as an 2d quantum field model on a ML which can be obtained by repeating 
the two rows in Fig. 2 in a vertical direction. This lattice is invariant under translations 
of the block of two R3 and four R2 matrices in horizontal direction as well as in time 
direction. Later we will show that the partition function of the model, Z = tr{T)^°, can 
be represented as a functional integral over scalar fermions as 

Z = tr(r)^« = J DiljDilje^-^'^+^^, (3) 

where the matrix Aij defines the hopping parameters along arrows on the corresponding 
ML. The matrix A inherits the translational invariance of the associated ML. The equiv- 
alence between the formulation in terms of scalar fermions as alluded to on the rhs of eq. 
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(jni) and the formulation in terms of the i?-matrices and r is demonstrated by passing to 
a coherent fermionic state basis and using the technique developed in [3 [TU] and will 
be given below. 

Let us recall the definition of the i?-matrices. The two-particle (i?2)jj-niatrix is an 
operator acting on the direct product of two two-dimensional spaces Vi and Vj with basis 
elements and |ji), respectively, as 

{R2U^l) ® = (-l)^(^^)(^(^'^)+^(^'^»(^.,)SiK2) ® IJ2) (4) 

and can be represented graphically as shown in Fig. la. 

Similarly, the three-particle (i?3)jjfc-matrix is acting on the direct product of three 
two-dimensional spaces Vi , Vj and with basis elements , and \ki) 

{R3hk\^l) ® ® |A:i) = {-l)(Pi^-)+P(nmPik^)+Pik.)H^^^^^ ^ |^-^) ^ |^^) (5) 

and can be represented graphically as in Fig. lb. In the expressions (@)) and (0) the parity 
factor p{a) takes into account the graded character of the elements \a) of the Fock space 
of scalar fermions 

n 

pia)=n, |«) = nc+|0). (6) 

i=l 

The basic states of the vector spaces Vj on each site i are \k), k = 0,1 with c|0) = 0, |1) = 
c+|0). 

In order to have a simple hopping model on the ML for both R2 and R3 we follow the 
construction of corresponding matrices in the SFM of the 3DIM on a DBBC lattice j3] 
and use as an ansatz an exponential of a quadratic form (see Q) of fermionic creation- 
annihilation operators: 

Ri=:expctAf^Cj : 1 = 2,3, (7) 

where the notion : (■) : means normal ordering for the odd i,j and anti- normal (hole) 
ordering for the even The matrix elements {Afj — 6ij), I = 2,3 coincide up to a sign 
with the hopping parameters Aij in Q 

Af^ - 6ij, i - odd, j - even, odd, . , 

Oij — Aij , I — even, j — even, odd. 



It is straightforward to calculate the matrix elements of Ri from ((Tj) in the Fock space 
basic. In particular we note that the fermion number is conserved: 
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The matrix elements [R^^^^^l^^ in the expression ()1U|) are connected with the Ajj's in 
(|7j) by the following equations: 
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The two-particle i?2-niatrix elements can be obtained from these expressions by taking 
A\f = A^j = everywhere. 

There are also some additional model dependent constraints for the R3 matrix para- 
meters coming from the specific lattice used: 
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These conditions express algebraically that we (by the lattice contruction) have no hop- 
pings across the hexagon (see Fig. 2), and this is the reason the i?3-matrix can not be 
reduced to the product of three -R2's, as one might have expected from general factoriza- 
tion properties. Going furter back to the original full SFM on a random ML coming from 
the 3DIM, the hexagon i?3-matrices are precisely associated with part of the embedded 
surfaces which carry a curvature, unlike the square R2 matrices which are associated with 
flat parts of the embedded surfaces. In the piecewise linear geometry it is well known that 
the curvature 

^n = ^(4-n) (14) 

is associated with the n-faces of two dimensional complexes. Therefore Kq = — vr/2, while 
K4 = 0. 

Viewing the R-matrices as associated with particle scattering one attaches, using the 
language of integrable systems, a spectral parameter to each of the particles, the spectral 
parameters being connected to the rapidities of the particles.. Therefore one expects that 
the i?3-matrix in general depends on three spectral parameters {p, r, s), while R2 depends 
on two spectral parameters (j9, r). 



3 Yang-Baxter equations and their solutions 

The integrability conditions for the model can be found in the standard way by defining 
and solving the associated YBEs equations. One can obtain the local YBEs using the 
graphical representation shown in Fig. 3a and Fig. 3b as described for instance in j3]: 

^i2(p, q)R234{p, r, s)Ri2{q, r)R34{q, s) = ^234(5, r, s)^23(p, r)R34{p, s)^34(p, q), (15) 

^i2(p, q)R23{p, r)R3A{p, s)Ri23{q, r, s) = Ruiq, r)Ru{q, s)i?234(p, r, s)'Ru{p, q), (16) 

These equations ensure that the transfer matrices T{p,r,s) = Tr{TiT2) (Ti and T2 are 
defined by Q) with different spectral parameters commute: 

[r(p,r,s),r(g,r,s)] = 0. (17) 

(we have just written the commutativity condition for one of spectral paramenters) . This 
is what we mean by the system being integrable since we can use fjl7|) to define an infinite 
set of mutually commuting, conserved charges Hn{r, s) by expanding the transfer matrix 
t{p, r, s) in powers of p: 

TT „N _ d^T{p,r,s) 

^n[r, s) — \p=p,. (18) 

We now look for non-trivial solutions to (|15p and (jl6|) where the intertwiner matrices 
Ri2{p,q), R^iip.q) have the structure ^ and (fT^ . For convenience let us use a specific 
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notation for the matrix elements of the R2- and i?2-niatrices, while keeping A\j s for the 
parameterization of -R3, ^3-matrices 

RooiP^r) =ai(p,r), R\l{p,r) =a2(p,r), 

Rli{p,r) =bi{p,r), R%{p,r) = -b2{p,r), 

RooiP, = ai{p, s), R\l{p, s) = aajp, s), 

Rllip, s) = hip, s), R%ip, s) = -h{p, s). 

With this notation the YBEs (fT3j) and (fTBj) reduce to the following constraints on the 

parameters a^, ak, bk, bk, k = 1,2. 

ai{p,r)a2{p,r) ^ ai{p,r)a2{p,r) ^ - 

b,{p,TMp,T) ^' h{p,r)b2{p,r) ' ^ ' 

ai{p,r)b2{p,s) = ai2{r,s), a2{p,r)b2{p, s) = a2i{r, s), 

ai{p,s)b2{p,r) = ai2{r,s), a2{p, s)bi{p,r) = a2i{r, s) 

and the R^ and -R3 elements are connected with a^, 6j by the relations 



^23 (9, ^^,5) = y{r,s), A^?(g,r,s) = x{r, s)ai{p, s), Af}{q,r,s) = v{r, s)a2{p,r) - 1, 
A23 {p,r,s) = y{r,s), A\{{p,r,s) = x{r, s)ai{p,r), A\^{p,r,s) = v{r, s)a2{p, s) - 1, 
Ai2{q,r,s) = w{r,s), Af^{q,r,s) = u{r, s)bi{p, s), A§^(g,r, s) = z{r, s)b2{p,r) - 1, 
A\2'{p,r,s) = w{r,s), Ay{p,r,s) = u{r, s)bi{p,r), Ay{p,r,s) = z{r, s)b2{p, s) - 1. 

(21) 

Here f{r),r{r), aij{r, s), aij{r, s),y{r, s),y{r, s),w{r, s),w{r, s),x{r, s),x{r, s), 
u{r, s),u{r, s), v{r, s),v{r-, s), z{r, s), z{r, s), are arbitrary functions of variables r, s. 



4 The Transfer Matrix 

In order to calculate the transfer matrix of our model explicitly, expressed as a normal 
ordered form of an exponential operator, one inserts ((Tj), given by (fT^ - (j7T|) . into 
(see Fig. 2). From Fig. 2 it follows that the double-row transfer matrix (0) is invariant 
under translations by four lattice spacings: we can restore the whole lattice structure by 
translation of the block of two R3 - and four i?2-niatrices, either horizontally or vertically. 
After some algebra, using Wick's contraction theorem, we obtain: 

t{p, r, s) = tr{TiT2) = F{r, s) : exp H{p, r,s) : . (22) 

The number-valued prefactor F{r, s) of the normal ordered exponential is 

F(r, .) = Af + - - , (23) 
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where N is the number of constituent horizontal blocks in the chain and 

t9 ± - 4£i£2 



Ai 



2 

= {1 - ai2io){l - a2iy) + {1 - uzai2a2i){l - vxai2a2i) - I, 

The operator H[j)^ r,s) in ()23|) is a quadratic form of fermionic creation and annihila- 
tion operators: 

H{p,r, s) = '^Hij{n — m)cf {n)cj{m). (24) 

This will be derived below. 

In general if is a nonlocal operator, Hij{n — m) being a polynomial function of Ai^2 
and £1^2 of degrees i — j, or N — {i — j), respectively. But by the translational invariance 
of the monodromy matrices we can Fourier transform the H[p, r, s) operators and present 
them in a compact form as 

N;A 4:;N 

H{p,r,s) = ^ Hij{n-m)ct{n)cj{m) = ^Hi,j{k)ct{k)cj{-k), 

n,m;i,j i,j',k 

1 ^ 

H^Ak) = T7Ee^'^^^..H- (25) 

n=l 

Let us now finally discuss how the calculation of the matrix elements of Hij{k) can 
be done by using the technique of fermionic coherent states with Grassmannian variables, 
as formulated in [SI |TT] . In this basis the product and the trace of the operators are 
obtained by integrating over the Grassmann variables. Choosing a anti-coherent and 
coherent basis for particles on the odd {2k — 1; n) and even {2k; n) boundary sites of the 
monodromy matrices (j21), respectively, (see Fig. 2), we define 

\^2k-i{n)) = e^2'=-i(")^i-i(")|0), (^2fc-i(n)| = {0\e^2,-i{n)^2,-,{n) ^ 

\i^2k{n)) = (c+ (n) - iJ2k{n))\Q), {i^2k{n)\ = {0\{c2k{n) + ^2k{n)) , (26) 

{i^2k-i{nM2k-i{n)) = e'^^fc-iWVa.-iW^ {^2k{n)\^IJ2k{n)) = e^^'^^'^)^^'^^"). 

These states are by construction eigenstates of the fermionic creation and annihilation 
operators and Ck with eigenvalues 4'k{n) and ipk{n) 

C2k I i'2k{n)) = -i)2k{n) I i)2k{n)) , {i'2k{n) \ = -{i'2k{n) \ ip2k{n), (27) 
ctk+i I i^2k+i{n)) = i'2k+i{n) I i'2k+i{n)) , {iJ2k+i{n) \ C2k+i = -{i^2k+i{n) \ iJ2k+i{n). 



7 



We attach also coherent states X2k{n), X2k{n), X2k+i{n), and X2k+i{n) to the inter- 
mediate sites between the two transfer matrices ti and T2. Then the full two row-transfer 
matrix in the coherent states basis, expressed via the one-row transfer matrices ti and T2, 
can be written as 

r{^,^) = J DxDxe~^^^Ti{i>2k-i{n),i/j2k{n),X2k{n),X2k-i{n)) (28) 
■^2(X2fc-i(n), X2kin),ip2kin),ip2k-iin)), 

where 

Ti = trTi, i = l,2 

N 2 _ 

= n n(^2fc-i(^)|(^2fe(^)|n|x2fc(^))lx2fc-i(^)), 

n=l k=l 

_ N 2 

T2(X,X,V^,^) = n '[[{X2k-i{n)\{x2k{n)\T2\ip2k{n))\ip2k-i{n)). (29) 

n=l k=l 

In order to define matrix multiplication in the coherent space basis we simply insert 
between operators the identity operators 

I dx^{n)dx^in)\x^in)){x^in)\e-^'^''^^^^''^ = 1, 

J dUn)dx^in)\Un)){x^in)\e-^^^''^^^^''^ = 1 (30) 

for coherent and anticoherent states, respectively. These relations simply express the 
completness of the coherent state basis. 

Now it is straightforward to express the transfer matrix in the coherent states basis. 
By inserting expressions (j3Up for the intermediate coherent states x between i?-operators 
in 121 and by considering matrix elements between external external quantum states ip we 
obtain 

t{iP,'iIj) = J DxDxe'^^^Y[Ri{n){'ilj,'ilj,XiX) = 
J DxDxDxDxexp{ (Xi('^)Aij(r2 - m)xj(m) + 'ipi{n)Aij{n - m)^/jj{m) 

n— m=0,l 

+i!i{n)Aij{n - m)xj{m) + Xi{n)Aij{n - m)tljj{m))}. (31) 

In this equation the matrix elements of the operator Ri{c,c) from ((7j) are represented as 

use of the coherent states. 
The matrix A represents vacuum fluctuations and its determinant will appear in the 
final expression after integration over internal states x- After Fourier transformation of 
the Grassmanian variables ip, x, the Fourier transformed A(/c) of the matrix A(n) becomes 



8 



a 10 X 10 matrix (10 is the number of the intermediate states in each repeating block, as 
one can see in Fig. 2): 
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The function F{p, r, s) in eq. can be expressed via the determinant of A as 
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while H{p, r, s) is defined by the Fourier transform of the matrix 

Hij{k) = Ai,{k) + A,Uk){A)-]{k)Kij{k) 
We will write down here Hij{k) only for a simple case: 
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and have found following matrix elements of the Hamiltonian 
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where 

A-^ = (1 - waf + (1 - u^a^) - 1 - 2 cos {2-Kk/N)u^wa^, (37) 

It should be mentioned that the dependence on the spectral parameter p in the expres- 
sions above comes from the fraction ai/hi. As seen we have a nonlocal model of hopping 
fermions. 

Finally we would like to make the following remark. In the article ^21 A. Zamolod- 
chikov has defined some S-matrix for scattering of straight strings and formulated the 
analog of the YBEs for them, called the Tetrahedron Equations, in order to have an in- 
tegrable model. In the vertex formulation this S-matrix ^3] has three initial and three 
final indices precisely as the i?3-matrix in our construction. Contrary to our situation, 
where particle number conservation is ensured by eq. ©, the non-trivial (to us known) 
solutions of the Tetrahedron Equations do not have particle number conservation, except 
of the one case represented in [Hj. 
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